Investigation of quantum transport by means of 0(N) real-space methods 
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We investigate quantum transport by a suitable development of the Kubo formula on a basis of 
orthogonal polynomials and using real space recursion approaches. We show that the method enable 
to treat systems close to a metal-insulator transition. On quantum Hall systems, results are given 
in the context of recent new universal relations between transport coefficients. RKKY mesoscopic 
interaction is also evaluated for two-dimensional quasiperiodic systems and insight of the usefulness 
of our method is given for large scale computational problems. 
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I. INTRODUCTION 



the associated eigenstates inferred a complicated nature 
described by multifractality. [pl]jl^ ] 
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Quantum transport in high dimensional non peri- 
odic systems, is usually investigated by way of scaling 
analysis through diagonalization of periodic hamiltonian. 
However, if N is the number of states, exact diagonaliza- 
tion requires a CPU time scaling as 0{N^), and memory 
scaling as 0{N'^). For sparse Hamiltonians, the use of 
Lanczos algorithm reduces memory and CPU time re- 
quirements to 0{N). jl]j^ In the present work, quantum 
transport in presence of magnetic field using a develop- 
ment of Kubo formula on orthogonal polynomials, is in- 
vestigated. This approach has been developed in the con- 
text of strongly disordered quasiperiodic Q and in- 
homogeneous systems where electronic susceptibility and 
conductivity were successfully estimated. [||j6j The key 
point of the algorithm is to rescale the density of states 
after an evaluation of the upper and lower bounds on 
energy and then to make a polynomial expansion of the 
associated Kubo formula. In principle, any orthogonal 
polynomials may be used, but it turns out that manip- 
ulations of Chebyshev polynomials are particularly suit- 
able as they are isomorphic to Fourier series. § Possible 
applications to nanostructures will be mentioned in this 
context. 



In two-dimensional systems, the study of localization 
in a magnetic field enables to address problems related 
to the quantum Hall effect and metal insulator transi- 
tions. Q A magnetic field in a pure system introduces 
further topological complication of the electronic spec- 
trum, which turns out to be a degenerate ensemble of dis- 
crete Landau levels. For disordered systems, in the limit 
of strong magnetic field, perturbational, numerical and 
field theoretical approaches have depicted a comprehen- 
sive view of the corresponding physical phenomena. |^] 
In real materials, due to disorder. Landau levels are en- 
larged and overlap to form Landau bands but extended 
states still persist at the center of each Landau bands, 
as revealed for instance by the divergence of the localiza- 
tion length near these critical energies. M The studies of 



The concept of multifractal states, also proposed by 
Kohmoto [T^ for quasiperiodic systems, is an impor- 
tant related issue. Indeed, anomalous quantum diffu- 
sion in quantum Hall systems [ |T^ (QHS) and in qua- 
sicrystals may have some direct relations with the 

observed physical properties. The possibility of an ex- 
perimental observation of multifractal exponents in the 
QHS has been proposed recently by Brandes et al. 
To understand properly the effects of these complicated 
eigenstates, one has to develop suitable methods. 

In this context, after discussing briefly the spectrum 
of the pure systems with magnetic field, we present a new 
method to evaluate the quantum transport coefficient in 
various situations. For QHS, the method brings comple- 
mentary results for studying the transition regions when 
Uxy goes from one quantized plateau to the next, while 
^xx goes through a peaked value and decrease again. It 
also enables a study of the phase diagram of the in- 
teger quantum Hall effect (QHE) through the calcula- 
tion of the transport coefficient as a function of disorder 
and magnetic strengthes. Besides, we show that physical 
phenomena indirectly related to electronic propagation, 
such as RKKY interaction, may be also investigated in 
high-dimcnsional aperiodic systems, by using similar al- 
gorithms. Finally, applications to large-scale computa- 
tional methods is suggested. 



II. NON INTERACTING ELECTRONS IN 2D 
SYSTEMS WITH MAGNETIC FIELD 

A. Spectral properties for pure systems 



The density of states on 2D disordered systems with 
magnetic field is a well established result which can also 
be investigated by recursion method, In the zero- 
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disorder limit, one gets a Landau level type spectrum 
where the number of gaps are defined by the dimension- 
less measure (magnetic strength) a — eBa^/2Trc (a: the 
lattice unit, B magnetic field). To investigate spectral 
and transport properties by recursion method, let's con- 
sider the tight-binding representation of the hamiltonian 



For n = 2, . . . , g and initial values ai = 1, /3i = 0, 71 = 
z,Si — —6^. The explicit form of the terminator and 
Green's function are thus given by : 



= ^ en\nx,ny){nx,ny\ + 

nx,ny 

A recursive construction of an orthogonal basis \^n) such 
that 

n\^n) = a„|*„) +6„|*„+i) +6„_i|*„_i) 

6„ = (*„+i|W|«'„) 

tridiagonalizes the hamiltonian and enables to evaluate 
the spectral properties quickly and accurately ^ . In the 
one-band case, the series a„,6„ usually converge after 
about 20 recursion steps, and the limit Ooc^oo give the 
final form of the Jacobi matrix, e.g. the convergence 
of the continuous fraction development of the Green's 
function is achieved. 
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Introducing a magnetic field with rational flux a = ^ 
leads to Q bands, separated by gaps. The asymptotic 
regime is given by q distincts values bq,q = 1 . . . , Q, 
moreover the terminator of the continuous fraction ex- 
pansion can be evaluated exactly : 



(V'olG(z)l^o) = 




z — ... 

so that the exact density of states can be determined. 
p8[ On Fig. 1, the recursion coefficients 6„ and the cor- 
responding TDoS for a = 0, 1/3, 1/8 are presented. The 
relation between the number of bands and the asymptotic 
behavior of recursion coefficients becomes obvious. 

We now introduce the disorder through the site en- 
ergies, which are chosen randomly (with uniform prob- 
ability) within the interval [— VF/2, M^/2]. The presence 
of disorder will smear out the gaps. We check that for 
W — 6t all the gaps have disappeared. 

In the following, after discussing current open prob- 
lems of the integer quantum Hall effect, the method we 
used to investigate Kubo formula ||l^, different from 
usual Landauer method pO|,pl| , is presented. It was ini- 
tially proposed by Mayou |3| for studying frequency de- 
pendent conductivity <t{uj) in disordered systems. An 
investigation of Kubo conductivity and anomalous quan- 
tum diffusion in 3D quasiperiodic systems has been car- 
ried out by developing a similar algorithm for static con- 
ductivity at zero temperature. B 



B. Effect of disorder on quantum transport in 
magnetic fields 



T(z) = ^ r(z) = "'^ + ^^^^^^ 

^ ' ^2 ^ ' ^ SqT{z) In the presence of disorder, extended states exist 

z only at energies for which localization length diverges as 

z i ^ ^ [E-Ec)-", V = 2.4±0.1. (9) At the critical value E^c, 

some controversies about the universality of the absolute 
z — bgT{z) value for the diagonal conductivity still remain. Experi- 

mentally ~ 0.2 — 0.5 which is in agreement with some 
^ith numerical scaling analysis of cF^^{E,rf). [p2|- 



7n — z^n-i bq^n+iOin-i^ oiyi — 7n-i Thc qucstion of the disappearance of these extended 

i/3„_i,/3„ = dn-1 states is important. In a study by Yang and Bhatt p4[, 
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a critical strength of randomness Wc by which all the 
extended states vanish, was found to be independent of 
magnetic field {Wc ~ Q.t, with t the constant hopping 
term). Since, in their calculations, Wc still persist in the 
weak magnetic field limit, the driven mechanism of the 
metal-insulator transition was consistent with the float- 
ing up picture of extended states, as proposed theoret- 
ically by Khmel'nitzkii and Laughlin. [^,^ According 
to their argument, the energies of extended states within 
a given energy range in the vicinity of the center of a 
landau band En = {n + l/2)hLJc tend to infinity 

if magnetic field is vanishingly small (as expected by An- 
derson theory of localization). This scenario turns out 
to be crucial for the global phase diagram proposed by 
Kivelson, Lee and Zhang (KLZ). ^ The KLZ approach 
leads to a physical relation between different quantum 
Hall liquids, known as the law of corresponding states. 
Accordingly, transitions from a quantum Hall states to 
insulators are allowed for certain values of = 1, 1/3, 
and forbidden for others = 2,3,4..., v = 2/5, ...). 

However, this picture has been recently contradicted 
by recent experiments by Song et al. where transition 
from ly — 2-state to insulator were observed. Sheng and 
Weng have also reported numerical support of a contin- 
uous disappearance of IQHE in a tight-binding model, 
and related to merging of extended states separating 
different plateaus. The authors claim that their re- 
sults in the weak magnetic field limit are not consistent 
with the floating up picture. In addition, they found 
that at the metal-insulator transition dissipative and 
non-dissipative conductivities are equal axx{Ec,Wc) = 
o'xyiEc, Wc). |p7| , p8| The physical explanation of this phe- 
nomenon still remains unclear and in the following we will 
show that their numerical calculations can be criticized. 



1. Kuho formula by recursion: diagonal conductivity 

The real-space calculation of the diagonal Kubo for- 
mula of the electronic conductivity may be considered 
as an alternative for usual Landauer conductance cal- 
culations, or diagonalization methods. The use of Lan- 
dauer formula for investigating quantum zero temper- 
ature transport is usually associated with free escape 
boundary conditions. One direction of the system is 
periodic whereas the other, of size L, is connected to 
metallic leads with different chemical potentials. Scaling 
analysis is performed through L. By recursion method, 
we avoid exact diagonalization of the Hamiltonian, so 
that we can treat in principle larger and more complex 
systems. To reduce the possible numerical instability 
at boundary conditions induced by the velocity opera- 
tor (periodic boundary conditions will indeed generate a 



short-circuit across the sample), we transform the Kubo 
formula in the following way {X{t) = f^-nt/n p^^-^m/n 
and X is the component along direction x of the position 
operator, the volume of the system) : 

\l t^oc t 

and we keep control of the asymptotic behavior of the 
quantum diffusion of the wave-packets. The conduc- 
tivity reads : 

^ I?,(t)x3m„^o < $,(i) I G{Ep+zr^) \ l,{t) > 

where | ^j{t) > = A'e"*^*/''| j > and | $j(t) > is nor- 
malized. The summation should be done over the total 
basis of states \jx,jy)i but it turns out that a limited 
number of initial sites is sufficient to achieve convergence 
of the calculation. The time-dependent evolution of a 
wave-packet initially localized at \jx,jy) is also evalu- 
ated by polynomial expansion of the evolution operator 
g-mt/n ^ dEVne-'^*^^)Vn{H) where we choose 

Chebyshev polynomials of first kind (see appendice). Q 

Finally, one could also define a scaling parameter 
from the construction of the recursion basis. Indeed, af- 
ter computing N recursion steps, the energy resolution is 
roughly F ~ ^ where W is the total bandwidth of the 
DoS. Starting from a localized state in I^^/q) = \jx,jy) 
the time corresponding to the "propagation" of the N 
states I^'tv) — VnI'^o) is r Accordingly, a scal- 

ing analysis as a function of N may be done. Besides, 
the use of Kubo formula implicitly requires taking the 
limit of an infinite system, as the spectrum of any fi- 
nite system is discrete. By retaining a finite ij imaginary 
part in Green's function, one replaces the delta function 
by a peaked smooth function of width 77, which must 
be greater than the level spacing in the finite system. 
The thermodynamic limit is achieved with increasing sys- 
tem size L — > 00, and 77 — > in order to retain all the 
contributions from the spectrum. In our calculation, we 
find that the thermodynamic limit is achieved for a finite 
number of initial states Jy). Concerning the transi- 
tion regions mentionned earlier, one notes that the effect 
of finite temperature or frequency is to smear the QH- 
metal-QH phase transitions. 

On Figure 2, the diagonal conductivities obtained 
for a disorder and magnetic strength W — 2,5 a — 1/3 
and r] = 0.15, are similar to those from reference. |Q The 
increasing of disordered strength leads to a reduction of 
the conductivity, as well as a shift of extended states to- 
wards the center of the band (as indicated by the arrows 
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on Figure 2). Ando first described numerically this ef- 
fect by studying the density of states and Thouless num- 
bers. [ p9|j30| A physical interpretation of this levitation 
mechanism associated to Landau-level mixing has been 
proposed for high magnetic field. 

The question of universality of the diagonal con- 
ductance at critical energies is also addressed in Fig- 
ure 3. As the finite imaginary part of Green's function 
tends to zero, one clearly sees that the e^/2/i-limit is ap- 
proached with our method. One notes that the figures 
(?7 = 0.15,0.09,0.05), in the inset, are obtained for only 
one initial state \ jx, jy), so that the fiuctuations of the ex- 
act shape are artificials. The central figure is an average 
result on ten different sites. 

We also discuss the numerical results obtained in HI] 
for a = 1/16 and different values of the disorder strength. 
To study the universal relation between the transport 
coefficient, the authors have computed (Jxx{Wc, Ep) and 
o'xyiWc, Ep) for Ep — —2.75 in proper units. On Fig. 

4, axx{W,Ep — —2.75) for different disorder strengthes 
is shown. The Fermi energy Ep = —2.75 considered 
in turns out to lie within a gap of the spectrum. The 
results from [ p3| are then inconsistent with the physics 
of the problem. Indeed, if the Fermi energy is located 
within a gap of the pure system, by increasing disorder, 
a consequent enhancement of the conductivity may be 
induced for a finite system, as it is illustrated on Figure 

5. However, if the Fermi level lies within the center of a 
Landau band, such a behavior is not seen in our results 
for the Fermi energy in the center of a Landau band. 
The authors of ||2^ may have not consider the correct 
location of the Fermi energy (remind that the position 
of the extended state will be affected by the disorder) 
and accordingly, their statement is at least numerically 
questionnable. 

On figure 6, we show the (Jxx{Ep) for two different 
values of disorder. The black curves are the results ob- 
tained for one site 'Dj{t) x Sm^^o < '^j(i) I G^Ep+irj) \ 
^j{t) >, whereas the red one is the averaged result. One 
can see that the fluctuations of the components, enter- 
ing in the sum of the Kubo expression in real space, is 
not critical, so that finite number of sites is sufficient to 
achieve convergence. 



to topological invariant known as the Chern numbers, 
which count the exact number of extended states up to 
the Fermi energy. Contribution to electronic conductiv- 
ity of a given eigenstate | k > can be calculated from 



ie'^h X - < k\Vy\q >< q\Vx\k > - < k\Vx\q >< q\Vy\k > 



' xy 



A 



■E 



< '^xy > 



dipid(p2(J^y{(pi,f2) = 



q - ^k) 
2 

A / 

fc 



On finite size systems, Hall conductance (boundary 
condition average) enables us to identify if < a'^y > is a 
so-called non-zero Chern number Mk or not, thus if the 
corresponoing |fe) is an extended state or localized state 
(formulation is due to Thouless, Kohmoto, Nightingale 
and den Nijs 's [ ^3|j3^ ). Recently, non-commutative ge- 
ometry has provided an interesting framework to inves- 
tigate Kubo formula and quantum Hall effect. |^ 

As we are interested, in particular, in computing the 
non-dissipative conductivity for energies in the transition 
regions, (where axy is not quantized), one has to adopt a 
different strategy. To that end, we propose a real-space 
approach of Hall conductivity. Starting from the general 
off-diagonal form of the Kubo conductivity, one shows 
that the proper algorithm allowing the expansion of the 
Hall conductance in a real-space basis is given by : 



2D. 



E2>Ef 



dEidE; 



fiE,)-fiE2) 
El — E2 



Ti[5{Ei^-H)y5{E2-n) 



enables us to expand the spectral measure on the basis of 
Chebyshev polynomials. After some simple algebra, one 
finds two parts to be evaluated separately : 



'xy 



^x: ^y/ 



where Imn is analytical and depends on the choice of 
the polynomial basis. In our case, it corresponds to 
(Ap =Arcos(-^^v^), and a.b associated to the weight 



function of Chebyshev polynomials) 



2. Hall Kubo conductivity 



The most spectacular result for 2D electrons in mag- 
netic field is given by the quantization of the Hall con- 

2 

ductance in units of when the Fermi level lies within 
a gap. This result has been well established and related 



^ 1 ( sm{m + n + 3)Ap sm{m + n + l)Ap^ 

Tib \ (to + n + 3) (to -|- n -|- 1) J 

The other part implies the calculation of the coefficients 
{jx,jy\'Pn{'H)yPm{'H)X\jx,jy). A reasonable number of 
initial sites \jx,jy) should be considered. The sum over 
m and n-indice is, given the form of the /„„ factors, lim- 
ited by some appropriate cut-off. One notes that the 
computational time is however much larger compared to 
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the Oxx calculation. Further work is in progress for test- 
ing the accuracy as well as the gain over conventional 
diagonalization procedures. |pa| 



III. RKKY MAGNETIC INTERACTION IN 
NON-PERIODIC SYSTEMS 



Amongst interesting phenomena related to elec- 
tronic propagation, the Rudermann Kittel Kasuya 
Yosida interaction (RKKY) between magnetic sites 
in disordered systems has been subjected to great atten- 
tion. In particular it has been shown to be very impor- 
tant for understanding the spin glass transition or more 
recently giant magnetoresistance effects in magnetic mul- 
tilayers. 

The RKKY interaction is generically given by 
'lRKKY{ri,rj,E) = j'^x{'ri,rj,E)Sri-Srj where J is the 
interaction between the localized moment S^^ and the 
spin of the itinerant electrons, and x(ri,rj, i?) contains 
the sum of all the electron-hole propagation paths from 
\ri) to \rj). The susceptibility can be written down as 



E' <E 



,{n\5{E~n)\r,){T,\5{E' 



E-E' 



Quasiperiodic systems cannot be described by 
such averaging process, so that the use of recursion 
method gives here some interesting quantitative informa- 
tions. §,§ On Fig. 7 and 8, the TDoS for a 2D quasiperi- 
odic Fibonacci quasilattice, as well as typical signature of 
aperiodic long range order are depicted. The strength of 
the quasiperiodic potential is Vqp — 0.4< (with t the con- 
stant hopping integral between first neighbors) and the 
susceptibility is given in a-units, with a the lattice spac- 
ing. In the Fig. 8, typical features of the interaction are 
shown, but due to the complicated nature of quasiperi- 
odic potential, no Fermi wavelength can be properly de- 
fined and oscillations exhibit resurgences that are absent 
from the periodic potential. 



IV. POSSIBLE APPLICATIONS FOR LARGE 
SCALE COMPUTATIONAL METHODS 



The expansion of any operator on orthogonal poly- 
nomials can also be applied in the context of large scale 
computational methods, which aim at reducing large 
memory and CPU time costs, for investigating more 
realistic models. For instance, spectral properties and 

ftic^l spectra for realistic model of Silicon quantum 
to na ve been evaluated and quantum confinement in- 
vestigated. [p8[ One can also consider thermodynami- 



and by development of spectral measure on Chebyshev 
polynomials one gets |^,^ 

X^J = 2 Ke^ Imnin I PniU) \ r,){r, \ P.^iU) \ r,) 

711, n 

_^,.P^{E)Pn{E') 



E>Ei 
E'<E, 



N{E)N{E')- 



E-E' 



dEdE' 



where the coefficients Imn have been previously defined 
for the Hall conductance. 

In metallic systems, the interaction is calculated ex- 
actly iRKKYir, E) ^ A{r) cos{2kpr + S{r))/r^, contrary 
to quasiperiodic or disorder systems where there is no 
simple analytical form. However, in weakly disordered 
systems, one can evaluate the quantum fiuctuations that 
arise in the higher moments of the interaction. jSTj It is 
found that only even moments lead to significant contri- 
butions : 



{x''{\n~r, I)) 



ix'^^Wn-r, I)) ^ exp( 



'P\ I 7-- _ rj 
n - rw 



with p{Ef) the DoS at Fermi level, the mean free path, 
whereas fip is constant independent of the parameter of 
the hamiltonian. 



cal properties of quantum systems thanks to the devel- 
opment of the partition functions Z{j3) = Tr[e~'^^] = 
5^„(/ dEe-P^Vn{E)) X VniU), with {Vn} a suitable ba- 
sis of orthogonal polynomials. Indeed, one needs to con- 
sider average quantities as {A) = Tr[pA\/Tr[p\ [p the 
density operator), but if the expansion of the partition 
function can be easily done in the case of a scalar argu- 
ment, when applied to an operator (for large scale sys- 
tems), one may not be able to switch into the eigen- 
representation of the Hamiltonian since full diagonaliza- 
tion is practically impossible. p9|-^ Consequently, poly- 
nomial expansions of operator are the only available al- 
ternative technique. 

Finally, in order to perform realistic atomistic calcu- 
lations of binding energies and interatomic forces, real- 
space methods turn out to be efficient and accurate 0(N)- 
methods. [53 



V. CONCLUSION 



Results on quantum transport by means of a real- 
space methods have been presented. The generality of 
the use of orthogonal polynomials for high-dimensional 
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non periodic systems has been discussed. In two- 
dimensional disordered systems with magnetic field, 
global phase diagram and new universalities of the quan- 
tum Hall effect may be studied by this method. One 
notes that, nowadays, lots of efforts are devoted to im- 
prove the efficiency of available numerical algorithms. 
The well-known Car-Parinello method has opened new 
perspectives for ab-initio electronic structure calcula- 
tions, but the development of the so-called "order-N 
scheme" |Q is of a major importance to investigate com- 
plex molecular systems. Quantum transport has also 
opened new challenges since the development of nanos- 
tructures (quantum dots, nanotubes,..) has become cru- 
cial in the field of electronic devices. In this context, 
further "order-N" schemes to investigate electronic prop- 
agation also need to be steadily improved. 



p{E)^ lim •^mG{E + iri) 

and general analytic properties of the Green's function 
can be investigated through the properties of the consid- 
ered polynomials. 
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VI. APPENDIX ON ORTHOGONAL 
POLYNOMIALS 



The principle of the real-space algorithms, used in 
this article, is based on expansion of spectral measures 
on orthogonal polynomials. Let us mention the require- 
ments as well as the key points of this approach. For any 
p{E) defined by p{E)dE — 1, there exists a family of 
orthogonal polynomials Vn [E) of degree n, such that 

/+^ N{E)Vn{E)r,n{E)dE = 6nm 

p{E)J2'Pn{E)Vn{E') = S{E - E') 

n 

for E, E' belonging to the spectral subset of p{E) and 
p{E) ^ p{E')- These polynomials generally satisfy 
a three-term recursive relation EPn{E) — a„7'„(_E) -|- 
hnVn+i{E) + bn-iVn-iiE) with 6_i = 0,n > 0, 
related to the moments of spectral function. How- 
ever, the expansion of spectral measure can be also 
done in a arbitrary basis of polynomials, for instance on 
Chebyshev polynomials Q„ : 

5{E -H)^ p{E) J2 Qn{E)Q„{n) 

71 

^ p{E)Y,'Pn{E)Pn{n) 



where p{E) = l/n^Ab'^ - {E - a)2 and p{E) the total 
density of states. Here, a and b define the band edges and 
bandwidth of the spectral function associated to Cheby- 
shev polynomials. Q Green's functions can be obtained 
by Hilbert transformation of p{E) (for absolute continu- 
ous spectrum) : 



G(z) 



dE 



P{E) 

{z-E) 
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VIII. FIGURES CAPTIONS 



Fig. 1. Recursion coefficients bn and corresponding to- 
tal density of states for three different strengtlis of tlie 
magnetic field. 

Fig. 2. Diagonal conductivity for quantum Hall systems 
with different disorder strengths {a = 1/3,7] = 0.15) in 
e^/h unit. 

Fig. 3. Diagonal conductivity versus Fermi energy and 
finite imaginary part of the Green's function r] (inset). 

Fig. 4. Diagonal conductivity versus Fermi energy and 
disorder strength for a = 1/16. 

Fig. 5. Dissipative conductivity as a function of disorder 
strength. Stars are associated with a Fermi level in the 

gap of the pure system (no disorder) whereas Fermi level 
lies within the middle of a Landau band for diamonds. 

Fig. 6. Components of the dissipative conductivity as 

a function of disorder strength and Fermi energy (black 
curves), and averaged result (red curve). 

Fig. 7. Total density of states for a two-dimensional 
quasiperiodic lattice {Vqp is the strength of the Fibonacci 
quasiperiodic potential). In the onset is given the energy 
interval for which corresponding susceptibility are calcu- 
lated. 

Fig. 8. Electronic susceptibility (in a- units, a lat- 
tice parameter) cLS cL function of Fermi energy with 
l.a,2.a,2.b,2.c,3.a,3.b,3.c corresponding respectively to 
Fermi energy -1.9 -|- 0.1 x A, A = 0, 1, 2, 3, 4, 5, 6. 
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